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Abstract — We consider the dependence of the internal structure of a 
neutron star in a close binary system on the semi-major axis of the binary 
orbit, focusing on the case when the Roche lobes of the components are 
nearly filled. We adopt a polytropic equation of state. The temporal 
evolution of the semi-major axis and its dependence on the mass ratio of 
the binary components and the polytropic index n are determined. The 
calculation are carried out right up to the moment of contact, when quasi- 
stationary model becomes invalid. We analyze differences in the shapes of 
the pulses of gravitational radiation emitted by a pair of point masses and 
by a binary neutron star, taking into account its internal structure and 
tidal deformations. 



1 INTRODUCTION 

Gravitational-wave astronomy is currently one of the most intensively developing 
branches of science. Coalescing neutron stars are especially interesting sources 
of gravitational radiation in the context of future gravitational-wave detectors. 
Binary neutron star exist in our Galaxy, and are observed as binary radio pulsars. 
Such systems lose energy via gravitational radiation; this causes the semi-major 
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axis of the system orbit to decrease, ultimately leading to the coalescence (merg- 
ing) of the stars. A large number of such systems for which the merging time is 
much less than the cosmological time are known, so that neutron stars merges 
should occur fairly often (no more rarely than about 3 events per year within a 
distance of 200 Mpc) [1-4]. Already operating and planned laser gravitational- 
wave detectors [3, 5, 6] will search for such events. 

The merging process can be divided into three physically distinct stages. 

1. The binary components move in quasi-Keplerian orbits without mass trans- 
fer, approaching each other only due to the emission of gravitational radi- 
ation. 

2. Upon the stars become sufficiently close, mass transfer between the binary 
components must be taken into account in addition to the Keplerian motion 
and gravitational radiation. 

3. In a late stage of the mass transfer, the merging of the binary components 
or falling of on star onto the other leads to the formation of a single ob- 
ject. This is the third stage of merging process. (Note that the distinction 
between the second and the third stages is quite arbitrary.) 

This problem has been treated in a large number of studies, among which 
two different approaches can be distinguished. The first approach is analytical. 
It is useful for studies of very distant pairs, in which distortions in the shape 
of the stars due to their mutual gravitational attraction and other tidal effects 
are negligible. In this case, the problem can be reduced to a form that permits 
finite, analytical solution. A number of such studies account for post-Newtonian 
gravitational effects [7-9] . This approach fails for even moderately close pairs, in 
which tidal effects, rotation, and the internal structure of the components become 
substantial, and must be taken into account when deriving information from the 
gravitational- wave pulses emitted [10]. 

The second approach relies on numerical computations. It uses hydrodynamic 
codes to describe the non-stationary coalescing and merging of two bodies, taking 
into account a whole series of physical processes that may be important. This 
approach is the only able to describe the late stages of a merger [11]. Unfortu- 
nately, current hydrodynamic codes have a number of drawbacks, including the 
presence of artificial viscosity in the numerical algorithms used and the inability 
to follow the evolution of a binary system from large values of the semi-major 
axis (say, ~ 10 6 stellar radii) up to contact. 

The first (analytical) method is able to trace the binary evolution up to a hun- 
dred stellar radii rather accurately, while the second (hydrodynamical) method 
can be used to calculate the evolution of the system beginning from several stellar 
radii. Thus, there remains an appreciable interval of semi-major axis values (say, 
from 100 to 5-10 R+) for which different methods must be applied. We describe 
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below one possible method in which quasi- stationary stellar configurations are 
calculated fully taking into account tidal effects. 

We will focus on a detailed description of the first stage, when the stars move 
only under the action of their mutual gravitational attraction and gravitational 
radiation. The problem will be solved numerically using the code developed 
in [12]. The second part of the problem is the determination of the shapes of 
gravitational- wave pulses emitted and identification of characteristic features that 
can be used to derive the main parameters of neutron stars. 

2 A MODEL FOR AN EQUILIBRIUM POLY- 
TROPIC STAR IN A BINARY SYSTEM 

We consider a model for a binary star with component masses M\ and M 2 and 
total mass M = Mi + M 2 . We assume the orbit to be circular with radius A; this 
assumption should be reasonable, since the orbit should be circularize long before 
the final stage of coalescence due to emission of gravitational radiation [7] . The 
primary component is treated as a polytropic star with polytropic index n, while 
the secondary is taken to be a point mass. Below, we will present arguments 
showing that this assumption is justified. As in [12], we chose the coordinate 
origin to be at the center of the first component, and the orthogonal coordinate 
axes to be oriented so that Ox passes through the centers of the two stars and 
Oy lies in the orbital plane. 

The equation of state of the first star is 

P = Kp 1+1/n 

where P and p are the pressure and density of the stellar matter and K is the 
polytropic constant. The density and the potential distributions are described 
by the equation of hydrostatic equilibrium: 

$1 + $ 2 + $ c + K(n + l)p 1/n = B , (1) 

where $1 is the potential of the first star (both inside and outside), $2 is the 
potential of the second star (point mass) 

GM 2 

\JJx - A) 2 + y 2 + z 2 ' 

$ c is the centrifugal potential, and B is a constant determined by the boundary 
conditions. The centrifugal potential is taken in the form 

<s> c = -y 2 n 2 [( x -x c ) 2 + y 2 }, 
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where X c is the center of mass of the system. In contrast to [12], we do not 
use the Keplerian formula for the semi-major axis for a given angular velocity Q. 
Instead, we determine this value self-consistently, taking into account the tidal 
interaction of the stars: 

^ (^i + ^ c )l ( ^o,o) = 0- 

On the other hand, the density and potential of the first star are related by 
the Poisson equation 

A$i = AnGp . (2) 

To find the equilibrium configuration for the extended component of the sys- 
tem, we must solve the system of equation (1) and (2) using the condition 
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pidVi = Mi . 



To implement the numerical algorithm, it is convenient to use the following 
dimensionless variables. We will express distances in units L = (GM/O 2 ) 1 ' 3 ) 
masses in units of M 2 , densities in units of M 2 /L 3 , the gravitational potential in 
units of GM 2 /L, and the polytropic constant in units of GM 2 1 ~ 1 / n L 3 / n ~ 1 . As a 
result, we arrive at the dimensionless equations 

0i + 2 + c + «(n + l)p 1/n = 6, (3) 

A0i = 47rp , (4) 

JpdV = q, q = M 2 /M 1} (5) 

a = A/L , (6) 



02 = - ; (7) 

(x — a) 2 + y 2 + z 2 



x c = q~ l J xpdV = , (8) 



c = - 1 / 2 (g + i)^^- TT ^J +y 2 j , (9) 

^(0i + 0c)| MiO ) = O, (10) 
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0i -> -q/\Jx 2 + y 2 + z 2 , x,y,z oo (11) 

in which lower-case letters denote dimensionless variables, k is the dimensionless 
polytropic constant, and p is now the dimensionless density. 

3 NUMERICAL ALGORITHM 

We transformed the system of equations to dimensionless form and solved it using 
the iterative method described in detail in [12]. We will now briefly summarize 
this method. 

Let Mi, M 2 , ft K, and n, be given, while A, B, p(x, y, z), and $i(x, y, z) must 
be determined. We do not assume a priori that Kepler's third law A 3 Q 2 = GM is 
satisfied. In dimensionless form, the input parameters are q, k, and n (note that 
rri2 — 1, mi — q, and uj = y/q + 1 in dimensionless form), while a, b, p(x,y,z), 
and <f>i(x,y,z) must be found. The dimensionless form of Kepler's third law is 
a — 1. All calculations were performed for the case q — 1 and n = 3 /2, and k was 
varied. This is in accordance with the relation K = kGM\~ x ^ u L^l n ~ x ^ which is 
equivalent to the shrinking of the system with fixed polytropic constant K, with 
k — > corresponding to an infinitely wide system (L — > oo). 

The calculations can be applied to systems with any M 1 and K. However, to 
make the results more transparent, we will scale them to correspond to a binary 
system in which M 1 = M 2 = 1.4M Q and the radius of unperturbed neutron star 
is 10 km. In this case of an infinitely wide pair, when the first component is 
purely spherical, this corresponds to an average density for the first component 
pi = 6.68 x 10 14 g/sm 3 and a polytropic constant K = 3.99 x 10 9 cm 4 g~ 2 / 3 s~ 2 . 

The scheme for the solution of the system of equation (3)— (11) is, in many 
ways, analogous to that in [12], with the addition of two more iteration cycles to 
solve the equations x c = and (10). These differences from [12] are associated 
with the necessity of determining a, which, in general, differs from unity (i.e. 
Kepler's third law is not necessarily valid for a distributed mass). We note, 
however, that the value for a returned by the computer code deviates from the 
unity by less then 0.5%. 

4 GRAVITATIONAL RADIATION FROM A 
BINARY SYSTEM WITH AN EXTENDED 
COMPONENT 

Energy losses due to emission of gravitational waves are described by the formula 
[13] 
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BE 
~dt 



45c 5 



G 




(12) 



where 



Qafs = P ( X * X P ~ l l^apx 2 ) dV 



Since the rate of change of the orbit at the end of its evolution is comparable 
to the velocities of the components, we must, in general, take into account the 
first and higher time derivatives of the semi-major axis. 

The main formulas for our analysis are presented for a rotating coordinate 
system. For transformation of quantities to a non-rotating frame, we write ex- 
pressions for quadrupole moments of the system and of individual bodies in non- 
rotating and rotating coordinate systems. We denote the non-rotating frame by 
Oxyz, and add a sign " over all quantities related to this frame. We keep the no- 
tation Oxyz for the rotating frame, and quantities in this frame will be presented 
without a sign " . The transformations from the rotating to the non-rotating 
coordinate system are given by equations 



where ip (t) = Q(t) is the time-dependent angular velocity of the binary system. 

Recall that our model for the binary system consists of a primary component 
that is a polytropic star under hydrostatic equilibrium in the rotating frame and 
a secondary component that is a point mass. The quadrupole tensor of such a 
binary system in the rotating frame can be written as 

Qaf3 = ■m x {r la r x p - 1 /z5 a { 3 r 2 ) + m 2 (r 2a r 2/3 - l / z 8 al3 rl) + q aP , 
where we have introduced the notation 



for the quadrupole tensor of the distributed mass, while the first and second 
terms are the contributions to the quadrupole tensor of two point masses located 
at r — r\ and r = r 2 respectively. Here, £ are the coordinates relative to the 
center of mass of the first star r — r\. 

It is known [13] that the quadrupole tensor can be expressed in terms of the 
inertia tensor of the first star 



x = x cos if) (t) — y sin ■?/>(£) , 



y = xsmip(t) + ycosip(t) , 





j >£pi(0(<Wf : - Up) 
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in accordance with the expressions 



where I = Ixx + Iyy + Izz is the trace of the inertia tensor. 

Since r is defined for the baricenter of the system, we can find explicit ex- 
pressions for the vectors of each component in the form 

m 2 A 

ri = j- 1, 0, , 

= _?£L(-1, 0, 0). 

Now, the explicit form of the quadrupole tensor in the rotating frame will 
be a function of two parameters: the quadrupole tensor of the system of point 
masses, which is determined by jiA 2 (// is the reduced mass of the binary system), 
and the quadrupole tensor of the extended component q a p. Since our model is 
symmetrical with respect to the planes y = and z — 0, only three of nine 
components of the quadrupole tensor are non-zero: qxx, Qyy, and qzz- Since the 
trace is equal to zero (qxx + Qyy + Qzz — 0), there remain only two independent 
parameters, which fully describe the total quadrupole tensor of the binary system 
in the rotating frame: 

_ fiA 2 q xx + qYY 
2 ' 

fiA 2 q xx ~ Qyy 

h ~— + 2 • 

We now write the components of the quadrupole tensor in the non-rotating 
frame: 

Qxx = /i + / 2 cos2^(t), 

Qxy = -/ 2 sin2^(t) , 
Qyy =fi- / 2 cos2^(t) , 
Qzz = -2/i , 
Qxz = Qyz = 0. 
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In the case of stationary orbit (Q— 0), the formula for gravitational radiation 
differs from the standard formula of two point masses [13] only in a term of the 
form q xx ~ Qyy- 

This formula reduces to the known relation [13] 

dE 32G , 2n2 32(7 fi „ 4 /MiM 2 



VP ; 5c 5 \ M J 
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5c 5 

in the case of point masses (gxx = Qyy = 0). Given the density distribution 
in the extended component of the system, we can calculate the dimensionless 
quadrupole momentum of this component with respect to the center of mass. 
Let us introduce the dimensionless function 



/ q 



l + q 



Qxx — Qyy 
M 9 L 2 



so that the luminosity of the binary system in the form of gravitational waves is 
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_d_E = 32GH4l^ 

Using the relation K = kGM^" 1 ^ X 3 /™ -1 , we eliminate L and obtain 

dE 32G^M^ _sn_ { q r , X „ x3 

-^T = 1 5 " « 3 -"h-T- + /( /c ) 1 + 9) 

dt 5c 5 R— \ 1 + 9 / 

Note that, in the case of binary systems in which extended components are 
elongated along the line through the stellar centers, taking into account the proper 
moment of inertia of a component will always lead to an increased energy-loss 
rate compared to the case of point masses. 

5 TOTAL ENERGY OF THE SYSTEM IN A 
MODEL WITH AN EXTENDED COMPO- 
NENT 

If the energy loss E is known, we can calculate the shrinking velocity of the stars 

E= dE/dA ■ A- In the standard approach [13], the energy of the system is taken 
to be the sum of gravitational energy of two point passes Ei 2 and their kinetic 
energies K\ and This yields an expression for the energy of the system 

E = —GMiM 2 /2A. If the masses are not point-like, we must also take into 
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account the self-gravitational energies E{ el ^ and E^ , the internal energies E\ nt 
and E™*, and the kinetic energies of rotation of the stars about their centers 
of mass K* and K 2 . We will take the second component to be point so that 
E s 2 elf = E l 2 nt = K* = 0. Then 

E = E 12 + K 1 + K 2 + Ef f + E[ nt + K{ , 
E 12 = J fft*dV x , E[ elf = V 2 J p^dVr , 

K{ = l / 2 Vt 2 I zz , E{ nt = nj PdV 1 = nK J p 1+1 / n dV 1 . 

Let us now consider the asymptotic of E[ nt , El e as L — > oo. In this case, the 
first component can be considered a polytropic, self-gravitating sphere. The so- 
lution of the equilibrium equations for such configurations are well known [14,15]. 
The total energy of a polytropic, self-gravitating sphere in hydrostatic equilibrium 
is given by the formula [14,16] 

.self ^ vint _ 1-nGMl 



E{ ai + Ef 



5 — n Ri 

The polytropic constant of such a sphere is determined by its mass and radius 
[14]: 



n— 1 3 — ri 



K = M n GM l " R x 

The coefficient M n is associated with the solution of Lane-Emden equationF] 
via the formula 

3 — n n — 1 

K = (4vr) 1 /™(n + l)- 1 e i 



W3/2 = 0.424 . 

On the other hand, 
lr The Lane-Emden equation has the form 

C 2 (fo')' = -o n , 

where £ is the independent variable and is the desired function (8' = d6/d£), with the 
boundary conditions given by 

6(0) = 1, 0'(O) = O, 
and by the values at the boundary of the sphere £1: = 0, and = — 
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We find as a result 



V 5 - n J 



E = G~M'^ n K'— . 
The sum of energies E 12 + K 1 + K 2 has the following asymptotic as L — > oo: 

E 12 + K 1 + K 2 = E (- l hq^) ■ 
We finally obtain for the total energy of the system the formula 

E = E (- 1 / 2 qK^ + h(K)) + . 

Here, h(n) is determined from calculations of this function (equal to the sum K\ 
and the differences of E[ eli + E{ nt and E 12 + Ki + K 2 from their limiting values), 
with h{n) — > as k — > (i.e., as L — > oo): 

, x E-EJ- l / 2 qn^) - 

h{K) = ^ * . 

Since the first component expands as the two components become closer, 
> 0. 

6 METRIC OF THE GRAVITATIONAL-RADI- 
ATION FIELD 

The physical observed quantities in gravitational-wave astronomy are the trans- 
versely traceless corrections to the metric tensor h 7 ^ and h 7T . These are ex- 
pressed in terms of time derivatives of the components of the quadrupole tensor 
of the system. To write the explicit form of the corrections to the metric tensor, 
we introduce a coordinate system that is rotated by the angle 9 about the Oy 
axis. We will calculate the transversely traceless corrections to the metric along 
the Ox axis in this new frame. The projection operator used to reduce the tensor 
to a transversely traceless for can be written 

/ \ 
P= 10 . 
\0 l) 

We can now obtain the transversely traceless components of the quadrupole ten- 
sor in the new coordinate system: 
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Ql T = 3 / 2 /i cos 2 9 - h 1 + f 26 cos 2^ , 

Q T J = f 2 sin sin 2^. 

The transversely traceless corrections to the metric are found using the time 
derivatives of the reduced components of the quadrupole tensor 

_ JG_ -"■ 



where Rq is the distance from the observer to the binary system. 

7 EQUATION FOR THE CONVERGENCE OF 
THE COMPONENTS; CALCULATION RE- 
SULTS 

The equation describing the convergence of a point mass and an extended com- 
ponent differs from that for two point masses. We will now determine how these 
equations change in our model and find the dependence of the main characteris- 
tics of the radiation on the parameters of the problem. 

First, some introductory words are necessary. The semi-major axis of the 
system is denoted A, while the main calculations use the arguments L. We 
emphasize that the difference between L and A is small. Therefore, we will 
write all functions as arguments of L, and assume A and L to be identical. The 
geometrical properties if the system for various k are depicted in Fig. 1, which 
shows the boundary of the star along the X-axis (asterisks) and the position of 
the point L\. The horizontal and vertical axes plot X and Li, both measured in 
terms of R 1 . In addition, the values if dimensionless polytropic constant k are 
marked on the axis to the right. We can see from Fig. 1 that the components are 
in contact if k = 0.153, which is corresponds to L/R 1 = 2.75 (these estimates are 
independent of the chosen values of M 1 and K, and are entirely determined by q 
and n). In addition, we can see from Fig. 1 that the size of the first component 
is very insensitive to the location of the second component when L/R 1 > 3.5: 

X m i n X max R\. 

Since the resulting formulas are quite cumbersome, we will restrict the follow- 
ing discussion to the case n = q = 1, for which we performed our calculations. 
We then have 
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4 . 5 i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i 
-2.0 -1.0 0.0 1.0 2.0 

X/R 1 

Figure 1: Geometrical characteristic of the system. The Ox axis shows the 
boundary of the extended component: the left mark corresponds to the far side 
with respect to L 1; and the right mark to the near side. R± is the radius of the 
extended component for an infinitely distant point-mass component. The solid 
line shows the position of the Lagrange point Li . The dimensionless polytropic 
constant k is indicated on the supplementary vertical axis. The contact of the 
extended component with L\ occurs at L = 2.7hR\. 



dE 256^M 2 2 °/ 3 2 
dt ~ 5cW (/i + lW) ■ 
A plot of the function /(/t) is given in Fig. 2, in which the solid line shows an 
approximation to this function using the formula f(n) = 171k 5A . We then have 
for the total energy of the system 



E _ G 2 M 7 2 /3 
K 



~ - 7 7 a/-3 /2 + • 



The function was also calculated. This function can be approximated by 
the relation h(n) = 13.5/t 4 - 1 (Fig. 3). 
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Figure 2: Plot of the function f(n) showing how much the intensities of 
gravitational-wave radiation differ for extended-component and point-mass mod- 
els. The solid line is the approximation f(n) ks 171k 5 ' 4 . 
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Figure 3: Plot of the function h(n) showing how the total energy of a system 
with an extended component changes as to the total energy for two point masses. 
The solid line is the approximation h(n) 13.5k 41 . 



13 







1 



3 



6 



10 15 20 25 30 



N, 



orb 



10 



-l 




2 



3 4 5 6 7 8 



9 



L/R 1 



Figure 4: Velocity of one of the component with respect to the center of mass 
of the system. The thick and thin lines correspond to the cases of extended 
and point masses, respectively The dashed line shows the orbital velocity of the 
system 1 / 2 AQ,. 

We performed the calculations for a system consisting of one point mass and 
one extended component (this was associated with limitations of our computa- 
tional capabilities, and this restriction will be relieved in the future). The results 
can be applied to two types of astrophysical objects: 

• binary systems consisting of a black hole and neutron star of the same 



• binary systems consisting of two identical neutron stars. 

The first case seems unrealistic because of evolutionary considerations, though 
there do exist models that predicts the formation of low-mass black holes in 
binaries [17]. 

We believe the second case to be more interesting. Using the above results, 
we can study such a system quite easily. Suppose the secondary component is 
also extended, and has the same configuration as the primary component, since 
their masses are equal. In this case the change associated with substituting the 
secondary point mass by an extended mass distribution is a correction of the next 
highest order, which can be neglected. It is necessary to double the corresponding 
terms in the formulas for the energy loss rate and the total energy of the system 
(f(n), h(K), and the term corresponding E^): 



masses; 
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Figure 5: Rate of radial shrinkage of the components r = Lj L in units of 
the orbital period. The thick and thin lines correspond to the case of extended 
and point masses, respectively. The dashed line shows the characteristic hydrody- 
namic time for the extended component, which indicates how quickly equilibrium 
is reached inside the star. The number of orbital periods until contact is indicated 
on the upper axis. 



dE _ 256G 9 M 2 20/3 5 



dt 5c 5 K 5 



-« 6 (y 2 + 2/(«)r , 



E = — j*- [-- - 2 3 / 7 A/" 3/2 + 2h(K) 
Differentiating the second formula, we obtain 

dE G 2 M 7 2 ,Z . nll . dn 

and finally 

dn 256G' 7 M 2 13/3 /t 5 ( 1 / 2 + 2/(/t)) 2 



dt 5c 5 K A l / 2 -2h'(K) 

Using the relations 

K 

GMo /3 L ' 
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dL GM 2 /3 L 2 dn 
~dt ~ K ~dt ' 
we can now obtain a differential equation for L: 




where r g = 2GM 2 /c 2 , which is r g = 4.15 km for the adopted value of M 2 . The 
auxiliary quantity a is defined by the relation: 

a(L) = 1 + 15.5 • (^i) 3 ' 1 + 13.3 • , (14) 

where R\ = 10 km is the radius of the unperturbed neutron star. Note that the 
maximum value a — 1.73 is reached at the moment of contact of the components, 
when L = 2.75R\. In the beginning of our calculations L = 8.5i?i, a = 1.02. 
The accuracy of this formula is better than 1% (this maximum value is reached 
at the point of contact). It is not difficult to verify that our results returns the 
standard expression 

dL 128G 3 Mf 
~dt ~ 5c 5 L 3 
for a point-mass binary when /(«) = h'(n) = and a — 1. 

Equation (13) describes the self-consistent evolution of a binary containing 
extended components due to emission of gravitational radiation. Our equations 
differ from those for a point-mass model in two new functions f(L) and h(L). 
Precisely these functions lead the differences in the evolution for these two models. 
The most important difference is the deviation of the function a(L) from unity 
in the model with extended components. This function is presented in a form 
convenient for comparison with post-Newtonian corrections to the law of motion 
for two point masses in a binary system. 

The evolution of a binary system with extended components differs most from 
the evolution for two point masses in the final stages, when the effect of tidal 
distortions becomes significant. Therefore, below, we will compare two cases: a 
binary system consisting of two extended components and binary consisting of 
two point masses. 

The next set of figures each show two lines: the thick line corresponds to the 
case of extended masses ant the thin line to the point-mass case. In addition, the 
number of orbital periods until contact is given on the upper axis on each plot. 

Figure 4 presents the right-hand side if differential equation (13) in the form 

of convergence velocity of the components V = x ji L (in units of velocity of 
light c) as a function of the distance between the components. The characteristic 
orbital velocity of the system 1 / 2 QL is shown by the dashed line. Note that the 
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Figure 6: Ratio of angular and centrifugal accelerations as a function of the 
distance between components for the cases of extended masses (thick line) and 
point masses (thin line). 

case of extended masses yields a higher velocity: at the end of convergence, this 
velocity is 0.016c for the case of extended masses and 0.005c for the point-mass 
case. Note also that the convergence velocity at the moment of contact would 
increase to 0.114c if we did not take into account both the mass distribution and 
the 'geometrical' distribution of the components; i.e., if we assumed that stars 
remain spheres with a fixed radius R\. 

Figure 5 presents the characteristic convergence time of the components r = 
Lj L in units of orbital period function of the distance between the 

components. The dashed line corresponds to the characteristic hydrodynamic 
time, which indicates the rate at which equilibrium is reached inside a star. 
We can see from this plot that the star as a whole is nearly in hydrodynamic 
equilibrium right up to the moment of contact. 

A comparison of the angular acceleration Q L and centrifugal acceleration 
Vt 2 L of the system is presented in Fig. 6 in the form of the ratio of these two 
quantities. This figure indicates the errors in our estimates of the power of the 
gravitational-wave radiation, which were obtained assuming that the components 
move along circular orbits at constant speed. 

Figures 7-9 depict the solution of the resulting differential equation for L in 
the form of the functions L(t), v gw {t) iy gw = 20/(271")) and Vg^ z {t). This last 
relation is interesting in that it should be linear in the case of point-mass binary, 
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Figure 7: Time dependence of the semi-major axis L(t) for model with extended 
masses (thick line) and point masses (thin line). 
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Figure 9: Linear gravitational-wave frequency to the — 8 /3 power for models with 
extended masses (thick line) and point masses (thin line). This quantity is a 
linear function of time in the point-mass case. 



since L~ L~ 3 ; L 4 ~ t — t; v gw ~ L~ 3 / 2 ). 

The waveforms h 7 ^ for the signals from two merging neutron stars are shown 
in Fig. 10. We can see that the stars in a binary with extended masses come into 
contact 10.5 ms earlier than in the point-mass model. The phase shift for these 
two cases at contact 5Af = 5.97. 

The time dependence of the correlation coefficient x between the signals from 
the point-mass and distributed binary systems is shown in Fig. 11. This correla- 
tion coefficient is determined by the formula 

\ 1/2 

X(t) = I h x h 2 dt I \ J h\dt J hjdt 
\o o 

where hi and h 2 are the signals from the extended-mass and point-mass systems, 
respectively. The correlation is disrupted 70 ms before contact. 

Figure 12 shows the binary system configurations (a) at the beginning of the 
calculation (t — 0, L — 85 km) and (b) at the time when the Roche lobes of the 
components are filled (contact). The time interval between these configurations 
id 17.8 ms, which corresponds to 32 revolutions of the system. 

Note that the Jeans criterion for the dynamic stability of the system [10] 



19 




Figure 10: Wave forms /i+ T for the signals from two merging neutron stars. We 
can see that a binary system with extended masses (thick line) reaches contact 
10.5 ms earlier than the point-mass system (thin line). The phase shift at the 
time of contact is 5Af = 5.97. 
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Figure 11: Time dependence of the correlation coefficients of the signals from 
extended-mass and point-mass systems. 
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Figure 12: Binary system configurations (a) at the beginning of the calculation 
(t = 0, L = 85 km) and (b) at the time when the Roche lobes of the components 
are filled (contact). The time interval between these configurations id 17.8 ms, 
which corresponds to 32 revolutions of the system. 
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J rot ^ / 3 J orb i 

remains undisrupted right up to the time of contact (J rot is the axial angular 
momentum of the components, and J orb is their orbital angular momentum). 

8 CONCLUSION 

After gravitational-wave detectors begin to operate, the first date on the fre- 
quency and characteristics of merging binary neutron stars will become available. 
The gravitational- wave signals from such events will be detected by convolving 
a known waveform (the calculated signal from two neutron stars in a point-mass 
model, for example) with the output of the gravitational- wave detector. The 
output of the gravitational-wave detector can be represented in the form 

s(t) = A ■ S(t) + n(t) , 

where S(t) is a given waveform with unknown amplitude and n(t) is detector 
noise. The search for a signal using a least-square method or other suitable 
technique reduces to finding the minimum of a functional with respect to the 
unknown parameter A, the signal amplitude. This procedure can be rewritten in 

terms of integrals of the form / S 2 (t) dt, and so forth. 

J o 

However, if the convolution is performed using a signal of the wrong form (in 
our case, the signal corresponding to two point masses) the effective signal will 
be diminished. This may result in false estimates of the signal level. Figure 11 
shows the convolution of the signal corresponding to our model with that for a 
pair of two point masses. We can see that the convolution with the standard 
waveform leads to the rapid decrease of the registered signal ~ 70 ms before 
contact. Note also that the region over which a(L) differs appreciably from unity 
falls into the frequency band for the peak sensitivity of planned gravitational- 
wave detectors [18]. Therefore, it is important to take into account the distortion 
of the waveform when attempting to register a gravitational- wave signal. 

It is evident that the distortions of the waveforms of gravitational- wave signals 
will impede detection of gravitational- wave events. On the other hand, precisely 
these distortions may prove to be very important for studies of the internal struc- 
ture and equation of state of neutron stars. 
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